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1 Introduction 



Recent years have seen remarkable progress in the theory of two dimensional quan- 
tum gravity, due in large extent to the discretized approach (see [Q, for a review). 
We have a good understanding of the situation for c < 1, where c is the central charge 
of the matter coupled to gravity, through the KPZ formula and associated re- 
sults. The situation changes for c > 1 where the KPZ formula gives meaningless 
predictions; numerical simulations f|] and analytical arguments || , |J , |7j favour the 
view that quasi one-dimensional configurations called branched polymers dominate 
the statistical ensemble in this region (there is, however, no definite proof of this 
conjecture). For generic branched polymer ensembles the critical exponent j s t r = \- 
Certain modifications of the weighting of different branched polymers lead to other 
positive values of 7 str (|§, ||). 

The KPZ results show that the critical behaviour of spin systems (e.g. Ising 
spins) interacting with two-dimensional quantum gravity is modified from that ob- 
served on a regular lattice. This raises the question of what would happen in an 
extended theory which interpolates between the fixed lattice (or crystalline) model 
and the gravity model ||10|| . The only two-dimensional problem of this type that has 
been solved is the case of R 2 gravity (here R is the scalar curvature density) |TT| . 
If surfaces of large R 2 are suppressed the ensemble becomes flatter and flatter at 
short distance scales. However it is found that 7 sir = — |, the pure gravity value, 
for all finite values of the R 2 coupling; only when all surfaces with non-minimal R? 
are completely forbidden does j str change. 

In this paper we will show that it is possible to change 7 s t r for a branched 
polymer (BP) ensemble by introducing matter fields, rather than fiddling directly 
with the geometry, just as matter fields at criticality can change j s t r in the two- 
dimensional models. Then we study the interpolation between "gravity" and "crys- 
talline" phases. 

As a first step, in Section 2 we classify the possible critical exponents obtainable 
from a generic class of quadratically coupled non-linear equations. In Section 3 we 
introduce and study the set of binary trees. Section 4 is concerned with the coupling 
of matter to binary trees. First we show that a normal Ising model coupled to binary 
trees yields the same results as pure trees. Introducing a more complicated model 
of Z 3 spins with a complex action leads to new behaviour, with j str taking values of 
2, | and | in addition to |. 
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2 Categorization of Critical Behaviour 



In what follows, we will be interested in the possible critical behaviour of systems 
of coupled equations of the form 

A i = z + (4 k A j A k . (1) 

The Ai are grand canonical ensemble partition functions and z is related to the 
chemical potential \i by z = e~ M . The coefficients of depend on the coupling 
constants of the specific model and the repeated indices are summed over; we assume 
that the set of equations is irreducible in the sense that no linear combination of 
AiS decouples. At small z the Aj are power series in z\ as z increases the Ai become 
non-analytic at some critical value z c . Defining 

Ma = 5 i5 - a\ k A k - at A, (2) 

and differentiating equation ( 15 ) w.r.t. z we obtain 

where w J = (1,...,1) and M _1 is the inverse of M, M _1 = det ^ M ^ (cof(M)) T . There 
are two ways in which the above expression can become non-analytic; either the Ai 
remain finite as z f z c but their derivatives blow up, or the Aj themselves blow up 
in that limit. We analyse of these two cases separately. 

If the Aj remain finite, but the ^ diverge, then we may write, with £ = z c — z, 

AiW^RiiO-C-^PiiC) (4) 

where 7 str < 1 and Ri(() and Pi(() are regular functions at £ = with Pj(0) ^ 0, 
Ri(0) ^ 0. Then 

^~r- (5) 

Now, since at 2 = z c the Aj remain finite, we must have det(M) = for the 
derivatives to diverge. The general expression for the determinant is obtained by 
substituting (f|) into (Q) and expanding which gives 

detM = C^^P^iC) + C 2 ( 1 -^)P( 2 )(C) + . . . + ^(l-7rtr)p(n)^^ (g) 

The -P^ n ^(C) depend on the a{ k and therefore on the coupling constants. 

We can now identify the different cases, starting with the most general; all the 
p( n \0) are nonzero. By comparing the leading powers in (Q) and (|5|) we see that 
j s tr obeys 

- Istr = "(I - Istr) (7) 
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and therefore 7 str = |. No constraints on the coupling constants are necessary for 
this case; so for example this value of j s t r can hold for areas of the phase diagram 
of a system with two coupling constants. 

The next case is that of P^(0) = but p( 2 )(0) 7^ 0. Then the expansion of the 
determinant (H) becomes 

detM ~ C 2 * 1 -***) + 0(C 2_7atr ) + 0{C 3(1 ^ atr) ) + ... (8) 

and the leading singularity of (H) is in the first term, leading to 

- Istr = -2(1 - 7 str ) (9) 

and to a value of 7 str = |. Notice that in order to make PW(0) = we are 
introducing a constraint in the coupling space. In the case of two coupling constants 
this leads to a line in the phase diagram. 

Further values of j str will be obtained if we introduce more constraints. If we 
demand P (1) (0) = P^(0) = but allow P^(0) ^ 0, the leading term will be 
the 0(C 3< - 1_7str - ) ) and we obtain 7 str = |; the two constraints necessary to achieve 
this situation mean that, in the case of two coupling constants, this value of 7 s t r 
is confined to a point in the phase diagram. A higher number of constraints will 
lead to other values of j s tr, but the above instances are enough for our purposes. 
Note that not all of the possible exponents found in this way can arise from any 
set of equations; the larger the number of coupled equations the more exotic is the 
possible behaviour. 

The second case is that the Ai themselves become infinite as z j z c ; of course 
the derivatives will also blow up. In this case 

Ai(z) = C^PiiO + RiiO (10) 

with 7 str > 1. Again detM must vanish at z — z c . Putting ( [TP] ) into (^) and 
expanding yields 

detM ~ C {1 ~ 7atr)mi C^ (11) 

i 

where m^rii are positive integers and a\ are (coupling constant dependent) constants 
such that 

- 7 str ) + rii > (12) 

so that we get 

- Istr = - Istr) + rij) (13) 

if a,\ = a 2 = ... = aj_i = and a,j 7^ 0. By assumption ^ str > 1 so that we must 



always have rii > 0. If = 1 then 7 sir = 1 (otherwise from (12) nii is negative 
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Table 1: Constraints and ^ s t r 



Istr 


Constraint 


Possible Values 


(n,m) 


# of constraints 


a 

2 


th - ?f > 


l 

2 


(1,1) 








1 


(2,2) 


2 






3 
2 


(2,1) 




5 
4 


^ - f > o 


3 
4 


(1,1) 








1 

2 


(1,2) 








1 

4 


(1,3) 


4 






1 


(2,4) 








5 
4 


(2,3) 




2 


rij — m.i > 


1 


(2,1) 


1 






2 


(3,1) 




3 


rij — 2mj > 


1 


(3,1) 








2 


(4,1) 


2 






3 


(5,1) 




4 


nj — 3rrii > 


1 


(4,1) 








2 


(5,1) 








3 


(6,1) 


3 






4 


(7,1) 





which is forbidden); this implies logarithmic singularities in which are not possible 
from a set of polynomial equations so we conclude that rij > 1. The question is now 
how many at must be zero (i.e., how many constraints there are) to get a certain 
value of 7 str . 

The number of constraints necessary to get different values of 7 str can be obtained 
by considering equation (|i~2|). For instance, if we demand j st r = §, (0) reduces to 
n i ~ f > 0. We can then examine the possibilites for n and m. Setting m — 1, we 
can put n = 1 (giving a value of |) and n = 2 (giving |); we can also put n = m = 2 
giving a value of 1. Thus to obtain j str — § we need a minimum of two constraints. 
The results of this analysis for some values of j str are shown in Table 1. 



3 Binary Trees 

The properties of the ensemble of tree graphs are well known |13[ . In this paper we 



will consider trees made of cubic vertices but modify the graphs slightly so that all 
the external lines except the root are attached to another line. This does not affect 
the recurrence relation generating the graphs which is shown in fig.|l|. Letting T/v be 
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Figure 1: The recurrence relation for tree graphs. 



the number of graphs with iV external vertices (not counting the root) we have 

JV-l 

T N =Y, T N-kT k , T x = 1 (14) 

k=l 



or 



OO 1 

T{z) = £ z N T N = - (l - VT^Ii 



(15) 



JV=1 



The exponent ^y s t r for the ensemble of graphs with one marked point (the root in 
this case) is defined so that the generating function G for the number of graphs of 
a given size has leading non-analytic behaviour 



dG 

dz 



1 Zrr Z ) 



~1str 



(16) 



so for the tree ensemble 7 str = |. Now we want to modify the problem so that 
one particular sort of graph is picked out and given a different weight; this is the 
ladder graph shown in fig.0. For a given iV > 3 there are precisely two of these, 




N vertices 



N=l 



N=2 




Figure 2: The ladder graph and the special cases for iV = 1,2. 

the one shown in fig.0 and its mirror image. However, there is only one graph for 
each of N — 1 and N = 2 (fig.0); we will define both of these to be ladders. All 
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the non-ladder graphs we will call "trees"; so in fact there are no trees for N < A. 
Letting be the number of ladders with N external vertices we have 



so that 



L x = L 2 = 1 
Ln>3 — 2 



and the exponent 7 s t r takes the value 2 for a pure ladder ensemble. The number of 



trees T' N satisfies 



T[ = T> 2 = T!, = 
T' N = q{ J2 ( T N-k + L N - k )(T' k + L k ) - L^-i - ^3} , iV > 3. 



'JV-l 



k fc=i ; 
The factor q enables us to assign a different relative weight in the ensemble to trees 
and ladders. A typical tree with iV external vertices gets a factor of q N ~ 1- , a tree 
which contains a ladder as a sub-graph has a smaller power of q associated with it. 
For the generating function we find 

T\z) = q{{T'{z) + L(z)) 2 - zL{z) - z 3 } (19) 

which is easily solved to yield the generating function for the modified ensemble 



G(z) = T'(z) + L(z) = - [1 - y/1 - Aq ((1 - qz)L - qz*)) . (20) 

For q = 1 this just becomes the usual tree generating function (p~9|) with j str = | 
whereas at q = it is equal to L(z). However for any positive non-zero value of q we 
find that 7 sir = |. This is easily seen by considering the behaviour of the argument 
of the square root as z is increased from zero; as z increases L(z) increases but 
before it diverges the argument of the square root must vanish (because it goes to 
—00 if L(z) goes to 00). Thus only at q = exactly do we manage to "freeze out" 
the general trees and get a system which contains only ladders. This behaviour is 
very reminiscent of the R 2 model discussed in the Introduction. 

Of course this phase structure is consistent with the general arguments we gave 
above. This is a system with one coupling constant, q, and for generic values it 
has 7 str = |; at one particular coupling, q — 0, it has j str = 2, which requires one 
constraint. 
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3 



Figure 3: Coupling of spins systems to a binary tree. A circle containing N repre- 
sents graphs made of trees with N external vertices excluding the root. 

4 Matter coupled to binary trees 

We can extend the model by coupling matter to the trees. Placing an Ising spin 
<jj = ±1 at each of the vertices we obtain the recurrence relation shown in fig.|3]. 
The interaction for two neighboring spins is given by e l3cri(Jj . Using the identity 
e P<Ti(Tj _ cosh^i _|_ tai(jj\ where t = tanh/3, we can take the link factor to be 
1 + tOiOy Then 

Y n-i 

^v = KlE £0- + t<Jia)(l + ta lC )(l + tbd)Z N ^ q (a,a 2 ,b)Z q ( Cl d,a z ). (21) 

q=l abed 

In the absence of a magnetic field the dependence of Zn on the external configuration 
a — (<Ti, cr 2 , a 3 ) must take the form 

Zn(ct) = A N + aia 2 B N + gig^Bn + g 2 <t 3 C n . (22) 

Inserting (^) into the partition function (^) and equating coefficients depending 
on the same combination of spins, we find a system of coupled equations for the 
An, Bn, Cn in terms of Atv-i, Atv-2, etc. Defining the grand-canonical partition 
functions 

oo 

A(z,t) = J2z N A N (t) (23) 

N=l 

and similarly for B(z,t), C(z,t), we obtain the system of equations: 

A = z + A 2 + t 3 B 2 (24) 
B = z + tAB + t 2 BC (25) 
C = z + tC 2 + t 2 B 2 . (26) 
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Singularities in any of the functions A, B or C (or in any of their derivatives) will 
signal a critical point in the full partition function (j2ip . 

The simplest situation is when the matter is decoupled, t = 0. Then B = C = z 
and A(z) is just the tree generating function T(z). The system can also be solved 
exactly for t — 1. Then the equations (§D-(§6|) imply that A — C = A 2 — C 2 which 
only has one sensible solution, A = C (the other, A = 1 — C does not exhibit the 
correct behavior as z — ► 0), which leads to A = B. Solving the remaining quadratic 
we find 

A = B = C= 1 -^^. (27) 
4 

Hence, z c — g and "y s t r = \- It is straightforward to show that 7 sir = ~ for all 
< t < 1. Eliminating A and C from (|2~^ -(|2~6|) gives 



B 



VI - 42 - 4t 3 F 2 + VI - 4zt - 4t 3 F 2 



f28) 



Clearly I? itself cannot diverge as 2 f z cr ; this is because the arguments of the square 
roots go negative when B is still finite. Hence j s t r < 1. Differentiating with respect 
to z we get 

fl£ B BH ( 1 t 1 
F— — = — + <^ = + } 29 

where 

2t 4 5 3 / 1 1 \ 

F = 1 = + = 30 

2 V VI - 4z - 4t 3 B 2 VI - 4zt - 4t 3 B 2 J v 7 

The right hand side of ( p9|) is clearly positive; the arguments of the square roots 
cannot vanish as z j 2 cr because the coefficient of ^ would vanish while they were 
still finite. So we are interested in the point where F vanishes. B is an increasing 
function of z and grows faster than z (by (|28"D ). Hence F is a monotonically de- 
creasing function of z at fixed t and there is a value of z for which F = and 
diverges. Furthermore |^ is a sum of negative definite terms; it follows that the 
generic value j str = \ applies everywhere and the cancelation necessary to get any 
other value of 7 sir cannot occur. 

The critical structure of the theory is not changed by the addition of Ising spins. 
This result is expected because the Ising spins never have a diverging correlation 
length in less than two dimensions and so cannot affect the global properties of the 
geometry. 

We now consider a generalised Ising model in which the spins take the values 1, 
e ±2 ^r WA - By allowing the partition function to include complex weights, we obtain 
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a richer phase structure than that of the 2-state case. The partition function is given 
by 

Z( a , P)=J2 e ^^)^^" 51 "^^^^)^^^*^ (31) 

s 

(note that in doing this it is necessary to assign an orientation to every link). The 
Boltzman factor for a link can be written in the form 

L = A(l + t{s\ Sj + a <a t))(l + 7 ( s t s . - Si s))) (32) 



where, 



A = 

7 = 



(l + 7 2 )(l-t 2 )(l + 2t) 
tan(i/3a) 



3 



< os(\/3q; 



* ~ 2cos( v / 3a) + e 3 ^' ( ^ 

We see from these expressions for that our model will be well defined for A in the 
interval [0, +oo]; however t must be restricted to [0, 1[ for A to be analytic. We can 
thus drop A and consider 

L = (1 + fiSiSj + uSjSj) (34) 

where /j, = t + 7(1 — t) and ^ = t — 7(1 — t) so that the partition function (|31~D 
becomes 

Z=^Z JV =^^n( 1 + + (35) 

trees trees {S^} links 

We can proceed as before to obtain a system of coupled equations which completely 
describe Z. The most general form of Z^, taking into account its global Z 3 symme- 
try, is 

Z N (Si, S 2 , S3) = A N + SiSlB N + SiSlB N + S\S2D N + S\S S D N + 

S2S3CN + S3S2EN + S1S2S3FN + S\sIsIGn. (36) 

By introducing this form into the partition function and then computing the sums 
leading to the grand canonical partition function, we arrive at the following set of 



equations: 



A = z + A 2 + /jiu(fjL + v)BD 

B = z + 11BA + fiuBC + u 2 fiDG 

B = z + fiBA + fi 2 BC + i/DG 
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C = z + fiCE + fivBD + v 2 fiFG 

D = z + uAD + ^DE + /j 2 uBF 

D = z + vAD + u 2 DE + f/BF 

E = z + isCE + fivBD + n 2 vFG 

F = z + ti{fi + v)EF + v 2 D 2 

G = z + u{n + u)CG + ix 2 B 2 . (37) 

For some curves in the (£,7) plane this set of equations can be solved exactly. We will 
find that, taken together with the general discussion of section 1, this is sufficient 
to enable us to elucidate the entire phase diagram. 

At t = 7 = 0, A is given as usual by the generating function of the Catalan 
numbers and all other quantities equal z so 7 str = \. At 7 = (this is the 3- 
state Potts model) we have that /U = v — t and the equations can be reduced by 
elementary manipulations to the set 

A = z + A 2 + 2t 3 B 2 

B = z + tBA + t 2 BC + t 3 BF 

C = z + tC 2 + t 2 B 2 + t 3 F 2 (38) 

F = z + 2t 2 CF + t 2 B 2 

B = B, D = D, E = C, G = F, 

which, like the Ising model (j24l) - (|26D , can be shown to have 7 s4r = | for all z. At 
t=l,(i = v= l independently of 7 and the system is again described by 
having ^ stT = \ for all 7; in this case the equations have elementary solutions 
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B = (39) 



16 4 



F = C 

and all the functions are non- analytic at z cr = ^. 

Along the line t = v = and the system of equations is greatly simplified 
and can be solved to give 

A = i(l_vT=4^ 
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B 



z(l - fiz) + [i 2 z 2 
1 — jiz — — [iz + fi 2 z)^ (1 — VI — 4z 



5 = 2 + 



1 — [IZ 

D = E = z 
D = z + f/BF 
F - 



1 — fl 2 Z 

G = z + ll 2 B 2 . (40) 



_ l 

str 2 
1 



str 



2 



At small fx the singularity is at z — \ and A, £>, S, G are non-analytic with 7, 
while C, D, E, F are analytic. So for small \i the u = line lies within the 7, 
region and is not a phase transition. However as /i increases we reach a point where 
the denominator of B vanishes for z < 4; this occurs at > /i where 

1 3 + d _ ^! = (41) 
4 P0 8 8 V ; 

which corresponds to /x — 1.2633, i.e. £ ~ 0.6316. So at /i = /x the singularity is 

at z — 7 but now 

5 ; 1 (42) 

and at this point j str = |. For fi > /xq we have that 

B — (43) 

with z cr < 7 so that A becomes analytic and 7 S f r = 2. From our general consider- 



ations we expect that j str = § occurs at a point and 7 str = 2 along a line in this 
two-dimensional coupling constant space. Hence for \i > /i the line v > separates 
two phases. 

It is straightforward to show from ([37])that for // 7^ 0, z/ 7^ 0, all the potentials 
are finite (this is done by assuming that a potential diverges at the critical point 
and obtaining a contradiction). Thus the two phases which are separated by the 
7 str = 2 line must both have j str = |, this being the only value which can exist on 
an area of the phase diagram of a system with two coupling constants. The nature 
of the phases is quite different. For £ > £0 and v > the thermodynamic limit exists 
in the usual sense, e.g. 

lim — log A N = const. (44) 

On the other hand when £ > £ and v < the potentials fluctuate in sign with 
increasing N; z cr has become complex. Of course this fluctuation is a remnant of 
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the complex (and sometimes negative) weights in the partition function. For t < to 
the line v = falls in the phase with the conventional thermodynamic limit as can 
be seen from the solutions ([TOp. The line separating the two phases must be in the 
v < region and, because all the potentials are finite, must therefore have 7 sir = | 
thus completing the phase diagram shown in fig.f|a. 

a) b) 

1 



1/2 3/2 


2\ 




2/3 - 




1/2 




Figure 4: The phase diagram for a) the complex Ising coupled to binary trees, and 
b) the complex Ising coupled to the ladder. Numbers indicate values of 7 str . 

Now consider the behaviour of the model when we weight ladder configurations 
differently from the rest of the ensemble as in section 3. At q = only ladder con- 
figurations survive; in the thermodynamic limit they are effectively one-dimensional 
objects, and it is straightforward to compute the transfer matrix to obtain the phase 
diagram depicted in fig. 4b. To elucidate the general structure of the phase diagram 
for the three coupling constants (u,fjL,q), one would in principle repeat the procedure 
that lead to equations (|I7| ) -(|2"0| ) . However it is clear that the general structure of 
the equations will be a system of coupled quadratic equations of the form discussed 
in section 2. Hence we can extend the analysis performed there to the present case. 

It is clear that the plane q = is exceptional; in this plane the coupled equations 
are in fact linear which is why j str takes integer values. Viewed as part of a three 
coupling constant {u^.q) space, y str = 2 can only exist on a plane and 7 str = 3 only 
on a line. On the other hand of the values for j str which occur at q — 1, 7 s t r = \ can 
exist on a volume, j str = | on a line, and •y str = | on a plane. If we move from q = 
to any positive q, we cannot have 7 s t r = 2 on a volume and hence the two regions 
which have that value at q = disappear for any positive q; the same happens to 
the line j s t r = 3. The general form of the phase diagram for q — 1 is mantained for 
any positive non-zero value of q. 

This situation mirrors the case of no matter coupling closely: there is an excep- 
tional plane q = where, by virtue of the "freezing" of the configurations, the values 
jstr takes are different from the ones at positive non-zero q. As soon as one moves 
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away from that plane, one is in a region where all configurations are important for 
the thermodynamics of the ensemble, and the results are qualitatively similar to 
those obtained at q = 1 where there is no damping of non-ladder configurations. 

5 Conclusions 

In this paper we investigated models of matter coupled to binary trees. Coupling 
the Ising model to such a system does not alter its critical properties, as might 
be expected but by introducing a complex-action Z3 model we showed that new 
behaviour can arise, giving values of ^ str different from ~. Interestingly, for a specific 
set of coupling constants v = 0, t > t , we find j str = 2 which is the behaviour 
expected when one single type of linear configuration dominates the full partition 
function; this is a sort of dimensional collapse driven by interaction with matter 
fields, although of a different sort from that studied in |J. On the other hand it is 
also possible to formulate these models to study explicitly the transition between a 
fixed configuration, the ladder, and the full fluctuating ensemble. We found that only 
when fluctuations are totally suppressed do the critical exponents change relative to 
the gravity case. 
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